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Abstract 


Turbularization  of  an  acoustic  boundary-layer  (Stokes  layer) 
impermeable  and  permeable  surfaces  is  analytically  considered. 


theoretical  approach- 


a  second-order  closure  model  of  turbulence. 


Both  an  approximate,  closed-form  solution  and  a  more  comprehensive  finite 
difference  solution  of  the  time  dependent,  parabolic,  one-dimensional 
governing  equations  are  obtained.  For  simple  acoustic  boundary-layers  on 
impermeable  surfaces,  both  the  approximate  solution  and  the  numerical  results 
for  the  critical  acoustic  Mach  number  required  for  turbulent  transition  are 
qualitatively  confirmed  by  experiment.  <~TheV  calculations  for  acoustic 
boundary-layers  with  transpiration  (injection)  indicate  a  substantial 
reduction  of  the  acoustic  Mach  number  required  for  transition,  up  to  a 
limiting  injection  velocity  that  is  frequency  dependent.  The  results  may 
provide  a  mechanism  for  flow-related  combustion  instability  in  practical 
systems,  particularly  solid  propellant  rockets,  since  turbularization  of  the 
near-surface  combustion  zone  could  result  at  relatively  low  acoustic  Mach 
numbers. 
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This  final  report  documents  a  completed  phase  of  work -performed  under  the 
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subject  grant,  the  scope  of' which /ye^nc€rns  the  analysis  of  turbulent  flow  and 
heat  transfer  behavior  in  rocket  chamber  flows  (C-systems).  Further  work  is 
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Foreword 


This  report  documents  a  completed  initial  phase  of  work  performed  under 
grant  AFOSR  85-0348,  the  subject  of  which  concerns  the  analysis  of  turbulent 
flow  and  heat  transfer  behavior  in  rocket  chamber  flows  (presently  termed  C- 
systems).  The  overall  scope  of  work  Includes  the  analysis  of  turbulence 
development  and  behavior  in  chemical  and  advanced  thermal  propulsion 
systems.  Specific  topics  addressed  include  combustion  instability  in  solid 
propulsion  systems  and  methods  of  heat  transfer  reduction  in  advanced  thermal 
propulsion  systems. 

The  initial  approach  to  these  problems  involves  assessment  of  their 
parametric  dependence  utilizing  a  well  characterized  second-order  turbulence 
closure  model.  Subsequently,  large-eddy  simulation  methods  will  be  developed 
for  these  highly  energetic  and  transitional  flows. 

This  work  is  continuing  under  a  new  grant  (AFOSR  86-0319). 
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1.  Introduction 


Oscillatory  flows  in  ducts  can  be  sustained  by  a  variety  of 
interactions.  These  range  from  purely  f luid-dynamically  induced  motions 
(e.g.  ,  large  scale  vortex  shedding  or  shockwave  instability),  to  the  more 
energetic  motions  possible  with  diabatic  flows.  Particularly  severe 
oscillations  can  occur  in  solid  propellant  rocket  chambers,  for  example, 
wherein  the  oscillatory  motion  can  be  driven  by  interactions  with  the 
substantial  energy  release  inherent  in  near-surface  combustion  processes. 

Prior  analytical  work  in  "combustion”  instability  in  solid  rockets  has 
identified  some  of  the  mechanisms  which  can  produce  velocity-field  coupling  to 
the  overall  instability  process  (see,  for  example,  the  review  by  Culick*). 
These  studies  have  shown  that  for  a  simple  longitudinal  standing  acoustic  wave 
in  a  duct,  a  Rayleigh  diabatic  instability  criterion  can  result  which  is 
dependent  upon  the  oscillatory  motion  of  the  gas  column.  Since  the  energy 
release  occurs  in  the  near  surface  region,  analytical  work  has  addressed 
acoustic  boundary-layer  effects  on  propellant  combustion  response  (see,  for 
example,  Lengele  ).  Implicit  in  several  response  function  analyses  is  the 
assumption  that  the  acoustic  boundary-layers  behave  quasi-steadily  and  in 
phase  with  the  longitudinal  acoustic  velocity  outside  the  boundary  layer. 

Recent  work  on  corabustion-f lowf ield  interactions  in  solid  rocket  chambers 
has  analytically  and  experimentally  examined  fundamental  fluid-dynamic  aspects 
of  mean  flow,  acoustic  wave  and  turbulence  behavior.  Hydrodynamically 
modeling  the  flow  with  a  semi-enclosed,  porous-walled  duct  (with  large 
injection  through  the  wall),  Brown  et  al^’^  experimentally  confirmed  the 
transitional  behavior  of  the  mean  flow  predicted  in  Reference  'S.  The  major 
emphasis  of  this  non-reactive  flow  experiment  was  to  investigate  the  effect  of 
low  and  moderate  amplitude  forced  acoustic  oscillations  on  the  flow,  using 
both  hot-wire  anemometry  and  surface-mounted  hot-film  sensors. 


Basing  their  conclusions  on  the  surface  sensor  measurements,  Brown,  et  al 
suggested  that  near  surface  turbulence  produced  by  the  mean-flow  transition 
process  appeared  to  "destroy"  the  coherent  response  of  the  sensors  downstream 
of  transition.  Consequently,  for  low  acoustic  amplitudes,  a  surface  response 
capable  of  inducing  instability  would  most  probably  originate  in  the  head  end 
(pre-transition)  section  of  a  rocket  chamber.  However,  it  is  noted  that  for 
large  axial  distances,  low  amplitude  acoustic  signals  appear  to  be  resurrected 
in  the  post  transition  region.  Lower  amplitude  harmonics  of  the  driving 
frequency  were  observed  by  Brown  in  both  the  surface  and  hot  wire 
measurements.  Their  existence  is  consistent  with  the  nonlinear  behavior  of  a 
(potentially  thick)  acoustic  boundary  layer  (Stokes  layer). 

Analyses  of  laminar  acoustic  boundary-layer  phenomena  involved  in 

propellant  response  have  been  offered  by  Flandro^;  Glick  and  Renie^;  Ben 
q  q  in 

Reuven  and  Hedge,  et  al  .  Baum  and  Levine  utilized  a  full  unsteady  Navier- 
Stokes  solution  method;  their  results  showed  a  Stokes  layer  persisting  even  in 
the  presence  of  strong  injection  velocities.  The  thickness  of  the  Stokes 
layer  was  of  the  same  order  as  the  Stokes  layer  for  noninjected  flow.  The 
presence  of  vortices  above  the  acoustic  boundary  layer  were  also  indicated  in 
the  results.  Although  several  aspects  of  these  analyses  are  interesting,  the 
following  important  point  is  noted.  Simple  scaling  estimates  of  the  laminar 
acoustic  boundary  layer  height  (the  Stokesian  thickness)  5  =/(u/f)  indicate 
that  is  well  above  the  gas-phase  flame  height  except  for  very  high 
frequencies  (f  >  10^  Hz).  Consequently,  minimal  acoustic  velocity 
interaction  with  combustion  would  be  expected  at  low  to  intermediate 
frequencies  as  indicated,  for  example,  by  the  more  detailed  analysis  of 
Reference  6.  Although  this  does  not  preclude  "velocity  coupled"  instability 
resulting  from  laminar  interactions  (whether  linear  or  nonlinear),  it  does 
suggest  that  other  mechanisms  of  interaction  should  be  explored. 


Toward  this  end,  the  potential  for  the  laminar  acoustic  boundary-layer  to 
undergo  transition  to  turbulence  is  considered.  It  is  well  established  that 
for  piston  driven  closed  duct  flows  (viz.,  in  which  there  is  negligible  mean 
axial  flow),  longitudinal  acoustic  waves  of  a  few  percent  relative  pressure 
amplitude  can  induce  turbulence  within  the  Stokes  layer.  The  literature 
review  and  independent  data  of  Merkli  and  Thomann^,  for  example,  confirm  a 
critical  amplitude  which  varies  as  /f  for  this  type  of  simple  acoustic 
motion.  Thus,  for  fixed  amplitude,  lower  frequency  disturbances  are  more 
likely  to  produce  turbulence. 

It  Is  known  from  studies  of  quasi-steady  propellant/ flowfield 

interactions  that  even  small  levels  of  turbulence  within  the  combustion  zone 

can  appreciably  enhance  propellant  combustion  rates.  Whether  acoustically 

induced  turbularization  can  occur  in  actual  rocket  chambers  depends  upon 

several  complex  effects  including  the  presence  of  large  injection  rates, 

surface  roughness,  and  of  course,  the  specific  type  of  combustion  process.  In 
1  9 

one  experiment  ,  values  of  the  threshold  velocity  for  propellant  response 
("acoustic  erosivity")  were  lowered  by  a  factor  of  2  relative  to  steady  state 
conditions. 

The  objective  of  the  present  analysis  is  to  consider,  under  several 
simplifying  assumptions,  the  general  trends  of  acoustic  boundary-layer 
turbularization  on  a  permeable  surface  in  the  presence  of  large  injection 
rates.  Though  pertinent  to  the  issues  discussed  previously,  the  posed  problem 
is  more  fundamental  and  of  potentially  broader  interest. 


The  motion  of  a  perfect  gas  described  by  the  Navier-Stokes  equations  is 
considered.  An  arbitrary  dependent  variable,  g(xt,t)  is  decomposed  according 
to  the  notation 

g(xlft)  =  g(xt,  t)  +  g'(xi,  t) 

with 


g  =  <g(xt ,  t)>  +  g'^Xj  ,  t)  . 


In  these  equations,  g'  is  the  turbulent  fluctuation,  g  is  the  ensemble 
average,  g"  is  the  acoustic  (deterministic)  component  and  <g>  is  the  long-time 
mean.  The  remainder  of  this  section  is  divided  into  two  parts,  the  first 
considering  order-of-magnitude  analysis  of  the  problem  and  derivation  of  an 
approximate  relation  for  critical  amplitude  of  transition  from  the  turbulence 
model  equations.  The  second  part  discusses  the  equations  and  method  used  for 
a  more  comprehensive  numerical  solution  of  the  acoustic  boundary-layer  in  the 
presence  of  transpiration. 

2.1  Order  of  Magnitude  Analysis 

Neglecting  third  order  correlations  and  the  axial  derivatives  of 
molecular  and  turbulent  stresses,  the  <  >  average  of  the  axial  momentum 
equation  is 

|-  (<p><u><u>+<p><u"u">+<u><p”u">+<u><p"u">) 


+  J—  ( <p><u><v>+<p><u"v">+<u><p’’v”>+<v><p"u">) 
dy 

3<p>  3_( <y>  5<u>)  SCO 

3x  3y  3y  3y 

where  the  mean  turbulent  shear  stress  is  <x>  =  -<p><u’v'>  (neglecting 
turbulent  density  correlations).  The  equation  for  the  acoustic  component  (u”) 


is,  to  first  order: 


<P>! +  (<P>u"  +  <u>p")||^  +  u"|h:<p><u>) 
(1)  (2)  (3) 


+  <p><u 


+  (<p>v" 


+ 


<V>p") 


3<u> 

3y 


(4) 


(5) 


u"37<<p><y>) 


+  <pXv> 


(6)  (7)  (8) 


*  &<■>$  +  £ 


(2) 


(9)  (10) 


where  r"  =  -<p>u'v*. 

In  equation  (l),  the  second-order  correlations  arising  from  the  unsteady 
motion  can  produce  virtual  stresses  and  convective  fluxes  analogous  to  those 
produced  by  the  turbulent  motion  alone.  These  correlations  induce  the 
classical  phenomenon  of  "acoustic  streaming,"^  wherein  a  component  of  the 
mean  flow  is  driven  by  acoustic  motion.  Of  potentially  more  importance  to 
propellant  combustion  are  the  analogous  second-order  terms  which  appear  in  the 
equation  for  the  transport  of  thermal  energy  in  the  propellant  flame  zone. 
While  these  thermal  transport  effects  will  not  be  analyzed  in  this  study, 
examples  of  their  importance  to  effective  steady-state  heat  transfer  may  be 


5 


found  in  the  literature  (see,  for  example,  Ref.  13).  The  effects  of  steady- 
state  enhancements  to  heat  transfer  from  acoustic  streaming  effects  have  not 
been  considered  with  respect  to  the  reactive  environment  in  rocket  engines, 
although  it  is  possible  that  these  effects  could  also  produce  a  coupling 
mechanism  for  instability. 

The  remainder  of  this  analysis  will  be  concerned  with  equation  (2)  and 

the  relative  importance  of  various  terms  within  the  acoustic  boundary  layer. 

To  further  simplify  this  equation,  it  will  be  assumed  that  the  density  in  the 

region  of  interest  is  approximately  constant  and  that  the  v"  produced  by 

propellant  response  is  negligible.  These  assumptions  would  necessarily  be 

removed  in  considering  a  comprehensive  propellant  velocity  response  analysis. 

The  magnitudes  of  the  terms  in  Equation  (2)  are  estimated  by  normalizing 

the  velocities  to  aQ,  the  chamber  length  to  L,  chamber  radius,  to  6  ,  and 

utilizing  the  continuity  equation  and  the  assumption  of  a  standing  first 

longitudinal  acoustic  mode.  It  Is  also  assumed  that  the  turbulence  shear 

2 

stress  may  be  scaled  by  x"  =  c^  <p>  u”  ,  where  1  c T  |  <<  1.  The  respective 

orders  of  the  terms  in  Equation  (2)  are  then 


{1,  <M> ,  <M> ,  <M> ,  (6  /5)<M>,  <M> ,  (a/f5  )<M  >,  1,  v/f6  a,  c  a/(fS  )  M"). 

Si  3  S  T  3  y  CUl  D 

The  order  unity  terms  1  and  8  recover  the  inviscid  acoustic  mode  solution 
in  the  central  region  of  the  duct.  A  considerable  simplification  results  if 
<M>  is  small,  since  terms  2  through  6  may  be  dropped  without  fundamentally 
altering  boundary  layer  characteristics  retained  in  terms  7,9,10.  (It  is 
Important  to  note  that  the  maximum  response  to  the  surface  sensors  in 
Reference  A  occurs  in  the  head  end  region,  where  <M>  is  <  0.1.)  Neglecting  for 
the  moment  the  transverse  convection  and  turbulence  terms,  it  is  seen  that  for 


the  viscous  term  to  be  of  order  unity  implies  the  Stokes  estimate, 


and  is  the  smallest  possible  laminar  boundary-layer  height  of  the  problem. 
The  effects  of  injection  on  the  acoustic  boundary  layer  (provided  largely  by 
term  7)  can  be  substantial  and  even  dominate.  The  ratio  between  the 
convection  and  viscous  terms  in  Equation  (3)  is  of  the  order  10  under 
conditions  of  interest. 

The  additional  shear  stress  provided  by  term  10  in  Equation  (2),  as  well 

as  the  analogous  additional  heat  flux  appearing  in  the  energy  equation,  result 

from  possible  turbulent  motions.  Further,  it  is  noted  that  term  10  is  the 

only  term  which  depends  on  the  amplitude  of  the  acoustic  motion.  (The 

acoustic  amplitude,  M" ,  is  a  linear  scaling  parameter  of  all  other  terms  but 

occurs  quadratically  in  the  turbulent  shear  term),  thus  posing  the  question  of 

a  possible  relation  to  nonlinear  stability  phenomena. 

As  noted  in  reference^,  there  are  a  few  approaches  which  can  be  employed 

to  obtain  a  stability  or  transition  criterion  tor  the  acoustic  boundary 

layer.  In  the  approach  adopted  here,  the  second-order  modeling  equations 

employed  in  reference  are  considered  (see  also  section  2.2).  Preliminary 

operations  and  assumptions  are  summarized  as  follows.  Equations  for  the 

velocity  correlations  u'^u'j  are  contracted  to  obtain  an  equation  for  the 

2  - 

turbulence  intensity  q  =  u'  u'  Fluid  material  properties  are  assumed 

1  X  • 

constant,  and  the  turbulence  (or  preturbulent  disturbance  level)  is  assumed 

small  so  that  third  and  higher  order  velocity  correlations  may  be  neglected. 

_  2  -  2 

It  is  also  assumed  as  an  approximation  that  u'u'  =  a  q  ,v'v'  =  a  q  , 

uu  vv 

2 

and  u'v'  =  a^q  >  where  aU(J,  avv  and  auv  are  constants.  Substituting  these 
relations  into  the  q^  equation,  taking  the  mean  flow  to  be  negligible,  and 
applying  boundary  layer  assumptions  for  a  quasi-planar  condition  yields: 


planar  condition  yields: 


A,3q  .  -  3q  .  Jq  ,  .  -  2.  3u"  3u" 

p(Tt  +  u  Tx  +  v  Ty  +  2p  q  (a  •  ~  +  a  -  ~ > 


uu  3  y  uv  3  y 


-2A  uq2  u32q 
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where  A  is  the  disturbance  raacrolength  -  scale  and  A  is  a  constant  used  in 
modeling  low  turbulent  Reynolds  number  dissipation.  The  approximation 
A  *  (where  Is  specified  by  Equation  (3)),  is  also  employed  to  evaluate 

the  maximal  effective  length  scale  appropriate  for  the  acoustic  boundary 
layer. 

•y 

Equation  (4),  which  is  linear  in  q  ,  may  be  order-of -magnitude  scaled  in 

2 

the  same  manner  as  Equation  (2).  The  requirement  that  Pj^j-  =  0  (for  neutral 
stability)  may  be  imposed  if  the  Reynolds  number  u"  6  / u  <  0(1).  The  order  of 
magnitude  analysis  indicates  that  the  fifth  (production)  term  on  the  left  hand 
side  and  the  first  (dissipation)  term  on  the  right  hand  side  of  equation  2  are 
dominant,  resulting  in  the  stability  criterion 
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Utilizing  the  turbulence  modeling  constants  specified  in  the  Appendix  of 
Reference  (3)  (i.e.,  A  =  3.25,  =  .17),  and  estimating  auv  =*  -.15  from  fully 
developed  flat-plate  turbulent  boundary  layer  flows,  yields  an  estimate  of  K  = 
750.  Merkli  and  Thomann  cite  prior  experimental  values  of  K  ranging  from  188 
to  915,  and  obtained  the  value  K  =  501  in  their  own  experiments.  (The 
constant  K  in  Equation  (5)  is  related  to  the  Merkli  and  Thomann  constant  Ac 


by  K  -  Ac/(2ir)/2)  .  They  speculated  that  the  variation  in  K  observed  in 

prior  studies  could  be  caused  by  variations  in  the  roughness  of  the  duct 
surface,  and  demonstrated  that  disturbances  caused  by  the  anemometer  probe  can 
also  appreciably  affect  K. 

The  critical  Mach  number  given  by  Equation  (5)  is  shown  in  Figure  1 
together  with  the  data  of  Merkli  and  Thomann.  As  can  be  seen,  the  functional 
dependence  is  correct,  and  even  the  very  approximate  estimates  of  physical 
constants  and  empirical  parameters  yield  quantitative  agreement  to  within 
several  percent  of  the  data.  To  the  authors'  knowledge,  the  derivation  of  an 
approximate  transition  relation  from  linearizing  this  type  of  complex 
turbulence  model  is  novel. 

2.2  Computational  Analysis 

The  second-order  turbulence  closure  approach  developed  by  Donaldson  and 
1  ft  1  7 

colleagues10’1  was  implemented  in  Reference  3  as  a  parabolized  model  for 
calculating  statistically  stationary,  compressible  transitional  flows  in 
porous  walled  ducts  with  large  infection  rates.  The  assumptions  and  order-of- 
magnitude  analysis  of  the  previous  section  again  yield  a  parabolic  equation 
system  for  the  present  acoustic  boundary-layer  problem.  Computationally,  the 
axial  convection  terms  pu  3g/3x  (where  g  is  an  arbitrary  dependent  variable) 
are  replaced  by  p  3g/3t  for  the  present  case.  With  the  exception  of  3p/3x  in 
the  momentum  equation,  all  other  axial  derivatives  are  taken  to  be  null.  The 
continuity  equation. 


3p  +  ]_  3_ 
3t  rv  3r 


(rV(pv  +  p'v' )] 
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and  the  3p/3t  term  on  the  right  hand  side  of  the  energy  equation  are  retained. 


however,  for  planned  work  involving  energetic  flows 


The  parabolic  differential  equation  system  may  be  considered  In  the 
functional  form 


at 
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where  g  =*{u,h,  u'u',  v'v',  w'w' ,  u’v',  h'u',  h'v',  h'h'}  ,  u,  v  and  w  are 

the  axial,  radial  and  circumferential  velocity  components,  v  Is  here  the  Index 
for  planar  or  axisymmetric  geometry,  p  the  density  and  h  *  cpT  Is  the  specific 
static  enthalpy.  The  molecular  transport  coefficient  p  represents  the  dynamic 
viscosity  (p)  or  the  thermal  conductivity  parameter,  a  *  k/cp,  as  appropriate 
for  each  equation.  However,  not  all  of  the  molecular  diffusion  terms  for  each 
equation  in  the  system  may  be  cast  in  the  form  shown  in  Equation  (2).  Those 
that  do  not  conform  are  implicitly  contained  within  the  complex  functions  Gg, 
which  also  represent  the  sources,  cross-coupling  and  dissipation  terms  for  the 
equations.  Although  the  more  general  (  )  averaging  Is  used  in  the  above 
notation,  the  only  mean  velocity  component  retained  is  v  ,  which  is  nearly 
equal  to  <v>  in  this  study  because  of  the  approximately  isothermal  and 
isobaric  conditions  assumed. 

The  pressure  along  the  duct,  p  =  <p>  +  p”(x,t),  is  specified  by  the  one¬ 
dimensional  standing  wave  solution 

p"  =  H  <p>cos(nwx/L)cos(nnf t)  (7) 

K13 

where  II  is  the  maximum  relative  acoustic  pressure  amplitude, 

ma 

and  u  =  II  a  /y.  All  calculations  to  be  presented  were  performed  at  the 

ma  ma  o 

velocity  antinode  of  the  first  longitudinal  mode.  Material  properties  were 
those  of  air  at  standard  conditions. 
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The  turbulence  length  scale,  A  ,  In  this  model  is  algebraic  and 
*1 

described  by  a  linear  variation  from  its  surface  value, 

Ag  (Ag=  0  for  smooth  walls),  to  a  plateau  level  proportionate  (with 
constant,  )  to  the  thickness  of  the  shear  flow.  For  the  present  analysis 
this  thickness  was  defined  as  6  ,  equal  to  the  height  above  the  surface  where 

the  boundary  asymptotes  to  99%  of  u"c.  For  injected  flows,  however,  waves  are 
convected  away  from  the  surface.  In  these  cases  the  first  zero-crossing  of  u" 
was  taken  as  the  effective  6  .  The  final  values  of  6  used  in  the  length 

dd  3d 

scale  expression  were  smoothed  by  integrating  in  time.  An  alternative 
differential  length-scale  equation  is  desirable,  but  much  has  been  written 
concerning  problems  with  the  low  turbulent  Reynolds  number  behavior  of  such 
equations. 

Boundary  conditions  at  the  duct  centerline  (or  centerplane)  are  the 
symmetry  conditions,  which  for  the  posed  system  are 

■I—  fu,  h,  u'u',  v'v’,  w'w* ,  h'h',  h'u'J  =  0 
o  r 

=  u'v'  =  h'v' 

Due  to  the  final  form  of  the  governing  equations,  the  requirement 
that  v  =  0  on  the  centerline  (centerplane)  is  necessarily  relaxed.  At  the 
duct  surface,  the  static  enthalpy  corresponding  to  a  given  temperature  is 
specified,  and  the  mean  injection  velocity,  vg  ,  is  prescribed.  The  no-slip 
condition  implies  that  all  tangential  velocity  components  and  their 
correlations  are  zero.  Experiments  have  confirmed  that  strong  injection  rates 
through  a  porous  plate  can  produce  pseudo-turbulent  disturbances  which  must  be 
included  in  the  ensemble  average  of  the  v’v'  boundary  condition.  The 
parameter  a ■  =  v'v' /v  was  introduced  in  Ref.  3  to  account  for  this 
important  source  of  disturbance.  All  enthalpy  correlations  are  null  due  to 


the  prescribed  uniform  surface  enthalpy. 

The  initial  profile  for  the  axial  velocity  component  was  taken  to  be  the 
analytical  solution  (at  t  =  0)  of  Sexl  and  Uchida  (in  Schlichting*®)  for 
oscillatory  flow  in  a  cylindrical  duct.  Initial  profiles  of  the  normal 
Reynolds  stresses  were  taken  to  be  isotropic,  and  proportionate  to  the  square 
of  the  initial  velocity  profile.  The  peak  initial  disturbance  level,  qmax(t  = 
0»y)/uma»  was  assumed  to  be  0.01  for  all  calculations  reported  here. 

The  implicit  numerical  procedure  and  adaptive  grid  are  described  in  prior 
studies.  Except  where  noted,  approximately  50  time  steps  per  period  were 
employed.  The  number  of  spatial  nodes  within  the  boundary  layer  varied  from 
about  40  for  laminar  flows  to  about  75  for  turbulent  cases. 


.  -V, >Vi  ‘  v  V 
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3.  Results  and  Discussion 


Normalized  profiles  of  the  axial  velocity  in  the  laminar  acoustic 
boundary  layer  (f  »  1000  Hz,  <S  =*  180  pm)  are  shown  for  various  times  within 
a  period  in  Figure  2.  The  results  display  the  classical  phase  shift  and 
Richardson  annular  effects  due  to  viscosity.  These  calculations  were 
performed  with  100  time  steps  per  period  and  produce  a  relatively  small 
computational  error,  principally  outside  the  boundary  layer  and  near  the 
maxima  in  the  acoustic  velocity.  Since  a  large  number  of  calculations  were 
needed  for  the  transition  studies,  the  number  of  time  steps  per  period  was 
halved  for  those  calculations,  and  the  error  in  the  calculated  acoustic 
velocity  was  then  approximately  doubled. 

Figure  3  shows  the  profiles  for  a  fully  (cyclically)  turbulent  boundary 

layer  at  100  Hz  with  II  »  .123.  There  is  a  noticeable  diffusion  of  the 

ma 

Richardson  effect  and  a  pronounced  increase  in  the  velocity  gradient  near  the 
surface.  However,  a  protracted  law-of-the-wall  (logarithmic)  region  is  not 
observed. 

Figure  4  shows  normalized  turbulence  intensity  profiles,  q/u^,  in  the 
boundary  layer  for  the  same  conditions  stated  for  Figure  3.  The  calculated 
peaks  in  turbulence  level  are  comparable  to  those  calculated  for  steady  state 
flows  (a  15%).  To  within  computational  accuracy,  these  peaks  are  symmetric 
with  respect  to  ±wft.  A  slight  phase  lag,  on  the  order  of  5-10  degrees,  exists 
between  the  intensity  maxima  and  the  rectified  acoustic  velocity  |u"c|.  This 
lag  Increases  at  higher  frequencies.  Also  note  that  the  intensity  decreases 
substantially,  hut  is  not  predicted  to  vanish  at  the  zero-crossings  of  u”c. 

Figure  5  shows  turbulent  velocity  profiles  for  a  Stokes  layer  with 
inlection  at  f  *  100  Hz,  II  =  0.080,  <v>  =  1  m/s.  A  pronounced  Richardson 


effect  is  evident,  surpassing  even  the  laminar  noninjected  results.  The 
estimated  value  of  the  effective  boundary  layer  thickness  for  this  case  was 
5.5  mm  (approximately  y/5fl  *  14  in  the  Figure).  The  convective  wave-train 
behavior  shown  was  also  evident  in  the  laminar  injected  Stokes  layer 
calculations  of  Reference  7. 

Figure  6  shows  the  corresponding  normalized  turbulence  intensity 
profiles.  Due,  in  part,  to  the  enhanced  Richardson  effect,  very  large 
intensity  levels  (3  302)  are  predicted.  As  turbulence  is  produced  in  the 
inner  layer  at  a  later  time  in  the  cycle,  the  velocity  wave  formed  from  an 
earlier  cycle  produces  a  harmonic  wave  in  the  intensity  profile  (ft  = 
5.125).  A  long  tail  is  predicted  due  to  large  decay  times  and  the  fact  that 
turbulence  is  beginning  to  accumulate  at  the  computational 
centerline  (5=5  cm)  where  the  symmetry  condition  is  employed.  In  an  actual 
duct  flow  environment,  the  mean  axial  flow  would  convect  these  tails 
downstream  as  they  interact  with  mean  flow  induced  turbulence. 

Calculated  results  for  the  critical  acoustic  Mach  number  at 
transition,  M"cr  =  H cr/Y »  are  shown  in  Figure  1  together  with  the  approximate 
relation  and  data  for  non-injected  Stokes  layers.  The  functional  dependence 
is  similar,  but  the  proximity  of  the  calculations  to  the  K  *  750  line  is 
coincidental.  Note,  for  example,  that  although  not  shown  here,  the  numerical 
results  display  a  dependence  on  the  initial  amplitude  of  turbulence  assumed. 
"Transition"  in  the  calculated  results  was  determined  by  monitoring  the  growth 
or  decay  of  initial  turbulence  from  the  initial  value  over  several  cycles. 
For  the  cases  with  injection,  the  growth  rates  were  not  as  strongly  dependent 
on  n  as  were  those  for  the  non-injected  cases. 

Critical  Mach  numbers  for  transition  as  a  function  of  the  normalized  mean 

V? 

injection  velocity,  <v>/(f\>)  z,  are  shown  in  Figure  7  for  three  frequencies, 


f  *  100,  300  and  1000  Hz.  A  pronounced  decrease  In  M“cr  Is  predicted  with 
increasing  injection  velocity  such  that  all  curves  exhibit  a  minimum  that  is 
frequency  dependent.  The  minima  are  almost  a  factor  of  three  below  the  M”cr 
for  noninjected  Stokes  layers.  As  an  example,  the  100  Hz  non-injected  Stokes 
layer  is  calculated  to  have  M"cr  ■  .083.  The  minimum  value  oi  M”cr  *  .025 


occurs  at  <v>  *  0.25  m/s,  while  M" 
of  the  line  at  this  frequency). 


.048  at  <v>  *  1.0  m/s  (the  termination 


Even  more  informative  are  the  same  results  for  transition  normalized  not 
as  a  critical  acoustic  Mach  number,  but  as  a  critical  acoustic  Reynolds 

lh 

number,  uma/(fv)  •  In  this  case,  all  three  curves  for  frequency  in  Figure  7 

collapse  onto  a  single  line  (to  within  computational  precision),  shown  as  the 

upper  line  in  Figure  8.  This  figure  also  shows  results  for  finite  levels  of 

the  surface  disturbance  parameter,  oy  =  0.035  and  0.07  for  two  different 

transition  criteria  to  be  discussed.  The  accompanying  value  of  ^  assumed  in 

s 

the  calculations  was  3x10”^  m.  These  values  are  viewed  as  realistic  and 
perhaps  conservative.  It  is  not  surprising  that  this  finite  disturbance 
effect,  which  continuously  "feeds”  the  injected  layer,  is  predicted  to  yield 
further  reduction  in  critical  acoustic  Reynolds  number,  analogously  to  quasi- 

*5 

steady  flows.  As  a  specific  example,  the  minimum  critical  Mach  number  for 

the  100  Hz  case  at  <v>  ■  .20  m/s  is  decreased  to  approximately  one  percent 

for  a  *  .035. 
v 

As  mentioned  earlier,  it  was  more  difficult  to  determine  the  value  of 
r  for  cases  with  injection  and  finite  levels  of  a  present.  The  first 
'terlon  employed  consisted  of  monitoring  the  maximum  turbulence  level, 
qmax/uma»  occur’-’ng  at  any  height  within  the  boundary  layer.  If  this  maximum 
attained  a  value  of  0.07  or  greater  within  3  cycles,  the  flow  was  considered 
turbulent.  Using  this  criterion,  the  ov  ■  0  lines  were  computed  along  with 
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the  alternating-dashed  lines  for  av  =  0.035  and  0.07.  While  this  criterion 

was  satisfactory  for  *  0  cases,  it  proved  to  be  unsatisfactory  for  cases 

with  finite  levels  of  o  . 

v 

To  illustrate  this  effect.  Figures  9a,  b.  ,  and  c  show  the  variation  of 
turbulence  intensity  versus  normalized  time  at  y  =  5  ,  for  three  values  of  II, 
<v>  =■  0.25  m/s  and  =  0.  It  is  evident  from  curves  9a  and  b  that  a  very 
small  increase  in  II  is  able  to  cause  the  flow  to  transition  to  a  cyclically 
turbulent  condition. 

When  finite  levels  of  a  are  included  the  determination  of  II  becomes 

v  cr 

more  complicated.  Figures  10a,  b,  and  c  are  results  for  <v>  =  0.25  m/s 

and  ov  =  0.035.  Figure  10a  is  a  case  for  which  the  turbulence  bursts  to  a 

q/Uma  level  of  approximately  0.24  but  then  dissipates  over  10  cycles  until 

only  the  low  level  turbulence  due  to  a is  still  present.  Figures  10b  and  c 

illustrate  the  increase  in  II  required  to  achieve  the  cyclic  behavior  which  we 

considered  to  be  characteristic  of  a  turbulent  case.  This  criterion  for 

transition  was  used  to  calculate  the  dashed  curves  for  a  =  0.035  and  0.07  in 

v 

Figure  8. 

To  recapitulate  the  results  of  the  two  transition  criteria,  the 
alternating  dashed  lines  in  Figure  8  qualitatively  represent  a  "burst"  or 
transient  criterion  for  transition,  while  the  dashed  lines  represent  the 
particular  (long-time)  criterion.  For  the  latter  case,  regions  to  the  right 
of  the  nearly  vertical  line  segments  are  indicated  to  be  turbulent.  The 

difference  between  the  two  criteria  is  in  part  explained  by  Equation  (4), 

3  2 

which  must  be  augmented  by  the  term  p  <v>  when  a ^  is  finite.  This  can 
result  in  three  different  turbulence  time  scales  being  competitive  in  specific 
parametric  cases 


16 


4.  Conclusions 


An  analysis  of  transitional  and  turbulent  acoustic  boundary  layers  In  the 
presence  of  strong  injection  has  been  presented.  The  problem  was  approached 
by  analyzing  the  behavior  of  a  second-ord  r  turbulence  model  rather  than  the 
traditional  Orr-Summerfeld  linearization.  The  approximate,  order-of-magnltude 
analysis  provides  a  simple  and  functionally  correct  estimate  of  transition  for 
simple  acoustic  boundary-layers.  This  technique  may  prove  useful  for 
estimating  the  stability  characteristics  of  other  types  of  flows. 

The  computed  results  for  Stokes  layers  with  injection  indicate  a 
substantial  increase  in  the  acoustic  boundary  layer  thickness  for  strongly 
injected  laminar  or  turbulent  flows,  but  no  "blow-off''  condition  is 
observed.  For  injected  fully  turbulent  flows,  a  pronounced  Richardson  effect 
is  obtained  and  is  accompanied  by  very  large  maximum  levels  of  turbulence  (> 
20%).  Both  axial  velocity  and  turbulence  profiles  exhibit  a  convected  wave 
train  shape.  The  turbulence  development  was  also  found  to  lag  the  rectified 
acoustic  velocity  by  a  few  degrees  at  lower  frequencies. 

Transition  resumes  for  conventional  Stokes  layers  are  in  agreement  with 
the  aDproximate  analysis  and  data  trends,  although  at  the  higher  acoustic 
amplitudes  shown  by  the  scaling  relation,  additional  nonlinear  effects  would 
be  important.  The  effect  of  injection  is  to  (non-monotonically)  decrease  the 
critical  acoustic  Mach  number  for  transition  by  up  to  a  factor  of  about 
three.  This  effect  is  frequency  dependent,  but  is  expressable  in  terms  of  the 
acoustic  and  injection  Reynolds  numbers  for  the  problem. 

A  further  appreciable  reduction  in  critical  Mach  number  (or  Reynolds 
number)  is  predicted  for  injection  velocities  with  finite,  continuous 
disturbance  levels.  This  effect  is  strong  enough  to  indicate  that  at  the 


minimum  critical  acoustic  Reynolds  number,  rather  modest  levels  of  acoustic 
pressure  ratio  (~  1%)  can  induce  significant  turbulence  levels  near  the 
surface. 


The  potential  for  turbularization  of  a  near-surface  reaction  zone  in 


ducted  flows  such  as  rocket  chambers  therefore  theoretically  exists.  However, 


whether  this  mechanism  can  produce  instability  in  such  systems  depends  upon 


several  effects  not  considered  in  this  investigation.  For  example,  combustic 


processes  and  (acoustically)  nonlinear  behavior  such  as  thermoacoustic 


streaming  could  influence  the  nature  of  the  overall  response,  while  the  axial 


mean  flow  could  bias  the  transitional  characteristics  over  the  acoustic 


velocity  and  pressure  cycles.  Further  research  on  these  topics  is  in 


progress. 
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7.  Nomenclature 


a  sonic  speed 

Cp  specific  heat  at  constant  pressure 

f  frequency,  Hz 

h  specific  sensible  enthalpy 

k  thermal  conductivity 

ks  equivalent  sand  roughness  height 

p  static  pressure 

i  '  Vo 

q  turbulence  intensity,  (u*  u^)  1 

qm  maximum  value  of  q 

R  inner  radius  of  a  cylindrical  duct 

Rec  axial-flow  Reynolds  Number,  ~P(^C  6/iIc 

Reg  Injection  Reynolds  number, 

Ret  turbulence  Reynolds  number,  ‘pTqA/u 

Ru  universal  gas  constant 

t  time 

T  static  temperature 

uj  velocity  vector  (u,v,w) 

Xj  coordinate  vector  (x,r,z) 

xQ  axial  distance  at  which  computational  Initial  conditions  are 

specified 

y  distance  from  surface 

distance  from  surface,  6-r,  (planar  flow);  or  R-r,  axisymmetric  flow 
characteristic  length  scale 
A  turbulence  macro-length  scale 

y  viscosity 
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V 


kinematic  viscosity 


density 
”  k/cp 

-  [v'v'/v2  |sf/2 
turbulence  shear  stress 
time  mean  of  variable 

Superscripts 

average  of  variable  over  turbulent  fluctuations 
*  turbulent  fluctuating  value  of  variable 

”  acoustic  fluctuating  value  of  variable 

Subscripts 

a  acoustic 

c  duct  centerline 

e  value  at  edge  of  boundary  layer 

h  condition  at  port  head  end 

m  max  absolute  value 

s  condition  at  surface 

,  differentiation 


P 

o 

av 

T 

<> 
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8.  Publications 


(of  work  supported  by  this  grant) 

During  the  period  of  performance,  the  following  publication  was  issued 
based  on  work  supported,  in  part,  by  the  subject  grant: 

Beddini,  R.  A.  and  Roberts,  T.A. ;  "Turbularization  of  an  Acoustic 
Boundary-Layer  on  a  Transpiring  Surface”;  Paper  No.  86-1448, 
presented  at  the  AIAA/ASME/SAE/ASEE  22nd  Joint  Propulsion  Conf., 
Huntsville,  AL,  June  1986,  (accepted  for  publication  in  the  A1AA 
Journal). 

9.  Personnel 

In  addition  to  the  Principal  Investigator,  Professor  Robert  A.  Beddini, 
the  following  graduate  students  contributed  technically  to  the  overall 
research  program  effort  and  were  supported  principally  or  secondarily  by  AFOSR 
funds  during  the  period  of  performance.  Degrees  being  sought  are  given  in 
parentheses. 

Principal  support  was  provided  to: 

Shyang-Lin  Kuo  (Ph.D.) 

Thomas  G.  Owano  (M.S.) 

Ted  A.  Roberts  (M.S.) 

Secondary  (summer)  support  was  provided  to: 


Claudia  M.  Bosch 


(M.S.) 


10.  Technical  Interactions 


During  the  period  of  performance,  two  principal  technical  interactions 
with  industry  representatives  occurred  and  are  presently  noted. 

1)  At  the  AFOSR/RPL  Rocket  Propulsion  research  meeting  in  September,  1986,  a 
technical  discussion  was  held  with  Dr.  Woodward  Waeche  of  Atlantic 
Research  Corp.  Dr.  Waeche  requested  the  Principal  Investigator  to  keep 
ARC  apprised  of  the  progress  of  the  work,  and  to  inform  him  should  any 
immediate  design  guidelines  for  control  of  solid  propellant  instability 
be  discovered. 

2)  At  the  AIAA  Propulsion  Conference  (Huntsville)  in  July,  1986,  a  technical 
discussion  was  held  with  Dr.  George  B.  Cox  of  Pratt  and  Whitney, 
Orlando.  Dr.  Cox  expressed  an  interested  in  the  work  since  P&W  utilizes 
porous  liners  for  passive  stability  control  on  several  of  its  liquid 
engine  combustion  chamber  designs.  He  requested  to  be  kept  informed  of 
future  work  that  included  effects  of  combustion,  heat  transfer  and  active 
cooling  of  the  chamber  surface. 
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Figure  1. 


Transition  to  turbulence  for  simple  acoustic  motion. 
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Figure  2 


Laminar  velocity  profiles  in  the  acoustic  boundary  layer 
without  Injection. 


Figure  3. 
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Turbulent  velocity  profiles  in  the  acoustic  boundary  layer 
without  injection. 
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Figure  4. 


Turbulence  intensity  profiles  in  the  acoustic  boundary  laver 
without  injection. 


>>■ 


1*1' 


n«o.o6 


1=100  Hz 
<v>»1  m/i 


VZ: 


f.f  =*.ooo 


Normalized  Distance  From  Surface,  y/<5 


Figure  5. 


Turbulent  velocity  profiles  In  the  acoustic  boundary  layer 


with  Injection. 
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Figure  6. 


Turbulence  Intensity  profiles  In  the  acoustic  boundary  layer 


with  Injection. 
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Figure  7. 


Figure  8. 
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Effect  of  Injection  on  the  stability  of  the  acoustic 
boundary  layer. 
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Effects  of  Injection  and  Injection  disturbance  level  on  the 
stability  of  the  acoustic  boundary  layer. 
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